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ABSTRACT 

Adopting the SPECTRUM package which is a stellar spectral synthesis program, we 
have synthesized a comprehensive set of 2,890 Near-InfraRed (NIR) synthetic spectra 
with a resolution and wavelength sampling similar to the SDSS and the forthcoming 
LAMOST spectra. During the synthesis, we have applied the 'New grids of ATLAS9 
Model Atmosphere' to provide a grid of local thermodynamic equilibrium (LTE) model 
atmospheres for effective temperatures (T e ff) ranging from 3,500 to 7,500 K, for surface 
gravities (log g) from 0.5 to 5.0 dex, for metallicities ([Fe/H]) from -4.0 to 0.5 dex, for 
abundances solar ([a/Fe] = 0.0 dex) and non-solar ([a/Fe] = +0.4 dex). So this synthetic 
stellar library is composed of 1,350 solor scaled abundance (SSA) and 1,530 non-solar 
scaled abundance (NSSA) spectra, grounding on which we have defined a new set of 
NIR Can triplet indices and an index CaT as the sum of the three. Then, these defined 
indices have been automatically measured on every spectrum of the synthetic stellar 
library and calibrated with the indices computed on the observational spectra from the 
INDO-U.S. stellar library. In order to check the effect of a-element enhancement on the 
so-defined Call indices, we have compared indices measured on the SSA spectra with 
those on the NSSA ones at the same terns of stellar parameters {T e ff, log g, [Fe/H]), 
and luckily, little influences of a-element enhancement has been found. Furthermore, 
comparisons of our synthetic indices with the observational ones from measurements on 
the INDO-U.S. stellar library, the SDSS-DR7, and SDSS-DR8 spectroscopic survey have 
been presented respectively for dwarfs and giants in specific. For dwarfs, our synthetic 
indices could well reproduce the behaviours of the observational indices versus stellar 
parameters, which verifies the validity of our index definitions for dwarfs. While for 
giants, the consistency between our synthetic indices and the observational ones does 
not appear that good. However, a new synthetic library of NIR Can indices has been 
founded for deeper studies on the NIR waveband of stellar spectra, and this library is 
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particularly appropriate for the SDSS and the forthcoming LAMOST stellar spectra. 
Later on, we have regressed the strength of the CaT index as a function of stellar 
parameters for both dwarfs and giants after a series of experimental investigations into 
relations of the indices with stellar parameters. For dwarfs, log g has little effect on the 
indices while [Fe/H] and T e ff play a role together and the leading factor is probably 
[Fe/H] which changes the stength of the indices by a positive trend. For giants, log g 
starts to influence on the strength of the indices by a negative trend for the metal-poor 
and even impact deeply for the metal-rich, and besides, [Fe/H] and T e ff still matter. 
Ultimately, a supplemental experimentation has been carried out to show that spectral 
noises do have effects on our set of NIR Can indices. However, the infection is weak 
enough to be ignored if the signal-to-noise ratio (SNR) falls below 20. 

Subject headings: galaxies:stellar content — stars:late-type — starsdundamental param- 
eters — techniques:spectroscopic — methods:data analysis — methods:statistical 



1. Introduction 

Due to more and more spectroscopic surveys, a huger number of spectra are turning available 
for studies on the physical properties (atmospheric parameters of stars; abundances, composition 
and ages of composite stellar populations of galaxies) by extracting information hidden in the 
absorption lines. However, problems of the limited sampling of spectral features caused by the 
detector pixel and the signal-to-noise ratio (SNR) seriously hamper the detailed fitting of the weak 
spectral features. So an alternative option is expected to extract physical information from the 
spectrum with no need of fitting faint absorption lines. Therefore, the concept of "spectral indices" 
has been introduced. This good idea for spectral analysis is to define a set of "spectral indices" each 
of which is an estimate of the equivalent width of a spectral feature measured by taking the ratio 
of the flux of feature in the central band to that of pseudo-continuum determined by two adjacent 



index systems of which the most successful is the Lick/IDS index system (Worthev et al. 


1994; 


Bur stein et al. 


1984; 


Faber et al. 


1985; 


Burstein et al. 


1986; 


Goreas et al. 


1993; 


Tracer et al. 


1998) 



which defines a set of 25 spectral indices in the optical waveband. Although it is widely used, 
Lick/IDS index system lacks for definitions of absorption features in the Near-InfraRed (NIR) 
waveband since the Lick spectra only cover from 4000 A to 6400 A in wavelength. 

In the NIR region (AA 7900 ~ 9100 A) of stellar spectra of types A ~ M, the strongest 
atomic feature is definitely Can triplet (AA 8498, 8542, 8662 A) which offers several advantages 
from the observational point of view for being the strongest metal line in the NIR spectrum, 
exceeding Nal (A8190) in strength by at least a factor of 5, and lying in a region re latively free 



of at mospheric absorptions. The Can triplet has been used as a luminosity indicator (IJones et al 



19841 ). and its strength is simultaneously a good metallicity indicator in stellar systems where the 
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NIR light is dominated by stars of a restricted range of surface gravities (jDiaz et al.lll989l ). Beside, 
it has provided a potential tool to study the properties of stel lar populations s ince the old stellar 
populations emit most of their light in this wavelength region ([Diaz et al.lll989l ). 



Progresses have been previ ously made at stud ies on the NIR Can triplet. By studying a sample 
of 62 stars of types FO ~ M5, IJones et al.1 (119841 ) has derived that the equivalent width (EW) of 
the NIR Can triplet only co rrelates stron g ly wi th log g in stars. This single- valued correlation 
has later been confirmed by ICarter et al.l (119861). However , by a set of definitions for the Can 
indices on a sample of 106 late-type stars, iDiaz et al.l (|1989i ) has demonstrated that in the high- 
metallicity range, the strength of the indices merely depends on log g while in the low-metallicity 
range, [Fe/H] turns to play a leading role. Afterwar ds, by refin e d defi nitions for the Can indices 
on a larger number of stars with a higher resolution. IZhou et al.l (Il99ll ) has clarified that T e ff has 
an important impact factor for th e low T e ff stars and log g plays a more important role for giants 
than dwarfs. ICenarro et al.l (|200ll ) has quantitatively estimated the strength of the Can indices by 
introducing a set of "generic indices". 



Although several works have been made on definitions for the NIR Can triplet indices (j Jones et al. 
19841 : biaz et al.l[l989l : ICenarro et aDboOll : Izhou et al.lll99lh . they are all based on samples of less 
stars with higher resolutio n and denser wav el ength sampling than t he Sloan Digital Sky Survey 
Data Release7 (SDSS DB7: Hfork et all fcoooh : kbazaiian et al.l (hood)') and the forthcoming Large 
Sky Area Multi-Oject Fiber Spectroscopic Telescope (LAMOST; iLuo et al.l (12008! )) spectra in the 
8300^4 ~ 8900A region. Besides, those stars also have inaccurate estimations for stellar parameters. 
In order to cope with these disadvantages, we are expecting to establish a new index system of 
the Can triplet which would be the most appropriate and convenient for use on the SDSS and the 
forthcoming LAMOST NIR spectra. We will perform this work on on a larger number of synthetic 
stellar spectra with approximately the same spectral resolution and wavelength sampling as the 
SDSS-DR7 and the forthcoming LAMOST spectra. 

The aim of this work is to constitute a new synthetic library of a set of the NIR Can indices 
for deeper studies on the NIR waveband of spectra, especially of the SDSS and the forthcoming 
LAMOST st ellar spectra. This lib rary could also be a supplemental part of the NIR indices to 
the work of IFranchini et al.l (|2010l ) which has established a synthetic library of several optical 
indices for SDSS and the forthcoming LAMOST spectra. In Section 2, we generate a large number 
of synthetic stellar spectra with various stellar parameters, on which a new set of the NIR Can 
indices has been defined and measured. Then, we also investigate the dependence of the Can indices 
on a-element enhancement. In Section 3, w e calibrate our syn thetic indices with the observational 
indices measured on the INDO-U.S. stars (IValdes et al.ll2004l ). Then, we compare the calibrated 
synthetic Call indices with the observational I NDO-U.S indices. Besides, tests of the calibrated 
synthetic indices on SDSS-DR7 and SDSS-DR8 (jAihara et al.ll201ll ) stars have been carried out. In 
Section 4, we analyse the relations of the calibrated synthetic Can indices with stellar parameters 
and approximately estimate the stength of the CaT index by stellar parameters. We compare the 
Call triplet indices with the CaHK index in Section 5. Finally, we summarize this paper in Section 
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6. As a supplement, we have performed a series of experiments in Appendix A to check into the 
effect of SNR on Call indices, and found that such an effect could be ignored until the SNR falls 
below 20. 



2. The Synthetic Indices 

2.1. Generate synthetic stellar spectra 

Above all, a saying is crucial to be stressed that it is not simple to model Can triplet lines using 
the local thermodynamic equilibrium (LTE) models since the line cores are formed in the stellar 
chromosphere and, he nce, non-local thermodynamic equilibrium (NLTE) models are required. So 



Jorgensen et al.l (|1992| ) has carried out, using NLTE models for the first time, the computation of 
the line cores, deriving full equivalent-widths and they have shown, however, that the effects of 
departures from LTE are negligible since the equivalent widths are dominated by the line wings. 
Considering that our Can triplet indices are estimated as equivalent widths, we could use LTE 
models for spectral synthesis without bringing departures in. Therefore, we compute the LTE 
model synthetic stellar spectra as our synthetic sample in this work, using the stellar spectral 



synthesis program SPECTRUM (jGrav et al.lll994l ). While using the SPECTRUM package, we set 



the wavelength sampling to be 0.05 (AA = 0.05), rotational velocity to be km/s, microturbulent 
velocity to be 2 km /s (e = 2 km/s), and pick the " New Grids of ATLAS9 Model Atmospheres" 



(jServen et al.ll2005l ) as the stellar atmosphere models from which we extract a subgrid as follows. 
Effective temperature(T e ^j): 3,500 ~ 7,500 K with a step of 250 K 
Surface gravity(log <?): 0.5 ~ 5.0 dex with a step of 0.5 dex 

For chemical composition, there are two assumptions which lead to different subgrids. One is 
the solar scaled abundances (SSA) grid ([a/Fe] = 0.0) which has 8 points for metallicity ([Fe/H] 
= -2.5, -2.0, -1.5, -1.0, -0.5, 0.0, 0.2, and 0.5 dex) and the other is the non-solar scaled abundances 
(NSSA) grid with [a/Fe] = 0.4 dex which has 9 points for metallicity ([Fe/H] = -4.0, -2.5, -2.0, 
-1.5, -1.0, -0.5, 0.0, 0.2, and 0.5 dex). After spectral synthesis of the above grids, we degrade the 
resolution of the synthetic spectra to R = 1800 (FWHM ~ 4.8 A in AA8300 ~ 8900 A) and change 
the wavelength spacing to be 2.0 A to match with SDSS and the forthcoming LAMOST spectra 
by using the SMOOTH2 task of the SPECTRUM package which convolves a Gaussian linespread 
function with the synthetic spectrum to smooth the spectrum to a desired resolution and then 
output the smoothed spectrum with a desired wavelength spacing. Hence, we finally derive 2890 
synthetic spectra composed of 1360 SSA model and 1530 NSSA model ones in the wavelength range 
of 8300 A 8900 A It is necessary to be emphasized that the shapes of the synthetic spectra of hotter 
stars (T e ff > 6000 K) appear a little abnormal for some reason so we have to choose T e ff = 6000 
K as a demarcation. We call the synthetic spectra with T e ff ^ 6000 K the Large Synthetic Sample 
(LSS) and those with 6000 K sC T eff sC 7500 K the Small Synthetic Sample (SSS). The LSS has 
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1870 spectra composed of 880 SSA model spectra and 990 NSSA model spectra while the SSS has 
1020 spectra consist of 480 SSA and 540 NSSA ones. 



2.2. Definition and measurement of the Can triplet indices 



An absorption strength is always expressed in terms of "an index " , which is usually composed 
of a "feature" bandpass and the flanking blue and red "continuum" bandpasses. The bandpasses 
should be defined by the following criteria: the "feature" bandpass needs to be defined to center 
on the feature of interest, the "continuum" bandpasses require to be located near the feature 
bandpass, to be in regions of less blended absorptions, and need for relative insensitivity to stellar 
yelocity dispersion b roadening. This last mandates a minimum length for bandpass definitions 
dWorthev et alJll994h . 



As is known, several atomic lines of intermediate strength (Fei AA8514.1, 8674.8, 8688.6, 
8824.21; Mgi A8806.8A and Til A8435.0A) exist in the wavelength region of the Can triplet lines. 
Moreover, the hydrogen Paschen series (AA8359.0, 8374.4, 8392.4, 8413.3, 8438.0, 8467.3, 8502.5, 
8545.4, 8598.4, 8665.0, 9750.5, 8862.8, from P19 to P10, respectively) are also blended into this 
wavelength region. However, the Can triplet absorptions are fortunatel y the strongest in s tellar 



spectra of types F5 ~ M2, and much stronger than other atomic lines (jCenarro et al.ll200ll ). Si- 



multaneously, the Paschen series nearly disappear in the spectra of stars colder than 6,0 00 K and 
thus could affect little on the Can triplet (see details in Figure 4 in lCenarro et al.l (|200ll )). 



To get a general idea, we check the centers, shapes, wings, and the local continuum of the 
Call triplet on the synthetic spectra of our LSS carefully. Then, we make definitions for the Can 
indices by choosing the central and continuum regions to cover the feature as complete as possible 
and simultaneously minimize the effect of the atomic or Paschen lines. By trial and error, we make 
a final definition for the Can triplet indices. The interval specifications for each line index are 
tabulated in Table 1. Like the Lick/IDS index system, each line index is measured by the ratio 
of the flux contained in a wavelength region centered on the feature relative to a pseudo local 
continuum represented by a straight line b etween the midpoi nts of the blue and red "continuum" 
regions close to the feature (for details see iFaber et al.l ()1985l )). Such a line-strength index whose 
definition is so similar to equivalent width gives a quantitative measure of a spectral signature. 
Additionally, we define a new compound index CaT as the sum of the three single index. 



Table 1: Wavelength intervals for each line index of the NIR Can triplet 



Identification 


Blue 


Line 


Red 




Bandpass(A) 


Bandpass(l) 


Bandpass (A) 


CaiiA8498^ 


8470.0~8490.0 


8491.0~8508.0 


8521.0~8531.0 


CaiiA8542i 


8521.0~8531.0 


8532.0~8553.0 


8563.0~8575.0 


CaiiA8662l 


8643.0~8653.0 


8655.0~8670.0 


8695.0~8710.0 
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According to the definitions, we automatically compute the Can indices on the LSS and SSS. 
These defined indices would as well be appropriate for studies of the NIR waveband of SDSS and the 
forthcoming LAMOST spectra. Conveniently, no consideration of matching resolutions is needed 
before using these indices on the SDSS and the forthcoming LAMOST spectra. 



2.3. Effect of a-element enhancement 



In this paper we investigate the effect of a-element enhancement on the Can indices. As 
our SSA model spectra do not have any a-element enhancement ([a/Fe] = 0.0) while the NSSA 
ones do ([a/Fe] = 0.4), we could check the effect directly via comparisons between the indices of 
the SSA spectra with those of the NSSA ones at the same terns of stellar parameters (T e ff, log g, 
[Fe/H]). For a clarification, we show the comparisons for each index in Figure 1 where it is obviously 
stated that the Call indices do not show any significant differences between SSA and NSSA indices. 
Furthermore, it seems that the consistency between the indices of SSA and NSSA is much better 
for the SSS (the red points) than the LSS (the black points); however, we can not be sure of the 
result for SSS because the spectra of SSS might not be correct from synthesis as is mentioned in 
Section 2.1. So we just believe that the NIR Can indices are all quasi- independent on a-element 
enhancement for stars cooler than 6000 K. 



3. Calibration of synthetic indices 



In order to calibrate our synthetic Can indices, we are required to compare them with those 
measured on reference stars with kno wn and accurate at mospheric parameters (T e fj, log g, [Fe/H]). 
Since the INDO-U.S. stellar library (jValdes et al.l 12004! ) covers the region of the NIR Call triplet 
in wavelength, we pick up stars with stellar parameters within our subgrid from the INDO-U.S. 
library as the test sample which includes 687 stars. Moreover, 564 (138 dwarfs and 426 giants) out 
of the 687 stars are cooler than 6000 K, so we select these 564 INDO-U.S. stars as the calibration 
sample. Beside s, the stellar parameters for the test sample we choose are from the version of 



Wu et al. 



(2011 



which has estimated pa rameters more a ccurately for the complete INDO-U.S. 
stars, particularly for giants (see details in lwu et al.l ( 2011 )). After degrading the resolution of the 
calibration sample to R = 1800, the 564 indices are computed on the resolution-matched sample. 
Here, we consider these 564 observational indices as the Observational Indices Sample (OIS). Then, 
we calculate the corresponding 564 synthetic indices with the same terns of stellar parameters as 
those of the calibration sample by linear interpolations in our SSA parameters- index grid (T e ff, 
[Fe/H], log g, index). Here, we regard these 564 synthetic indices as the Synthetic Indices Sample 
(SIS). By performing linear regressions between OIS and SIS, we derive a set of transformation 
coefficients (slopes and constants). Expectedly, the slopes of the linear fits are all close to one, which 
gives us confidence on the consistency between theoretical and observational spectra. However, the 
constant factors for dwarfs differs from those for giants. In detail, the constant calibration factors 
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are respectively 0.196, 0.645 and 0.425 for Cal, Ca2, and Ca3 index for dwarfs and 0.340, 0.719 
and 0.555 for giants. As an example of calibration, we show the calibrations of the three indices 
for dwarfs in Figure 2 where the black lines represent the calibration lines which could transform 
our synthetic indices into observational ones. Similarly, the transformations for Cal, Ca2, and 
Ca3 index are respectively Y = 1.238X + 0.340, Y = 1.258X + 0.719, and Y = 1.148X + 0.555 
for giants. As the Can indices are all quasi-independent on a-element enhancement (See details 
in Section 2.3), we do not repeat the calibration process for NSSA synthetic indices. Instead, the 
derived sets of transformation coefficients are used to convert both the SSA and NSSA synthetic 
indices into the observational indices system. 



3.1. Comparison between calibrated synthetic and observed INDO-U.S Indices 

In this section, we compare the calibrated synthetic indices with the observational indices 
measured on the test sample from the INDO-U.S. library. We use T e ff as the b aseline for com 



paris ons since T e ff is the most important and less model-dependent parameter (jFranchini et al 



20101 ). The temperature scales of the synthetic and the observational spectra are different which 
might introduce mismatches between the two scales, so we need to assume a temperature scale 6 
= 5040 /T e ff for both of the synthetic and the observational spectra, and we check the validity 
of this assumption by studying the behavior of the synthetic and the observational Call indices. 
Additionally, we divide the sample into two groups in terms of log g. One group is "dwarfs" with 
log g ^ 3.5 and the other group with log g < 3.5 is "giants". Then, we perform the comparisons 
in different [Fe/H] bins seperately for the groups. 

For each of the groups, the values of the Call indices computed on the test sample from the 
INDO-U.S. library are plotted versus 9 in eight different [Fe/H] bins which have intervals of A [Fe/H] 
= ±0.1 dex respectively centered at -2.5, -2.0, -1.5, -1.0, -0.5, 0.0, +0.2, and +0.5 dex. Then, in 
each bin, both SSA and NSSA isogravities of the calibrated synthetic indices are overplotted for the 
sake of comparisons. It should be noted that we revise the interval to be A [Fe/H] = ±0.3 and ±0.2 
when [Fe/H] centers at -2.5 and 0.5 dex for the group of dwarfs since there would be no suitable 
stars in these two bins if the interval still maintained to be A [Fe/H] = ±0.1 dex. Wheares, for the 
same reason, the interval is changed to be A [Fe/H] = ±0.2 when [Fe/H] centers at -2.5 -1.0, and 
0.5 dex for the group of giants. 

Then for dwarfs, Figures 3, 4, 5, and 6 show the spectral index Cal, Ca2, Ca3, and CaT versus 
9 (6 = 5040/T e ff) respectively in the eight different metallicity bins. In each of the bins, synthetic 
SSA isogravities are presented by solid lines and NSSA isogravities are shown as dashed lines while 
the observational indices from the INDO-U.S. library are presented by black dots. We draw a black 
dashed line vertically at T e ff = 6000 K in every bin to distinguish the indices of hot stars (hotter 
than 6000 K in the left region to the dashed line; Hot Indices Region) from the indices of cool stars 
(cooler than 6000 K in the right region to the dashed line; Cool Indices Region). In Hot Indices 
Region, the observational point positions might not match well with the synthetic isogravities 
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Co1 Ca2 Ca3 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 1 2 3 4 5 6 1 2 3 4 5 

SSA SSA SSA 



Fig. 1. — The effect of a-element enhancement on the Can triplet indices. The x-axis presents 
index values of the SSA spectra ([a/Fe] = 0.0) while the y-axis shows those of the NSSA spectra 
([a/Fe] = +0.4). The black and red points stand for the LSS and SSS spectra, respectively. 




Fig. 2. — Calibration of synthetic NIR Can triplet indices for dwarfs. This plot shows values of 
observed indices versus those of the corresponding synthetic indices at the proper terns of atmo- 
spheric parameters. The black lines with gradiants close to 1 represent the transformations which 
could transform our synthetic indices into the observed indices system. The black lines for Cal, 
Ca2, and Ca3 is respectively expressed as Y = 1.061X + 0.196, Y = 1.000X + 0.645, and Y = 
0.932X + 0.425. 
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somewhere because of the possible irregularities of the synthetic spectra of SSS which indicates 
that the synthetic isogravities in Hot Indices Region are perhaps not what they really are. In Cool 
Indices Region, it is obviously shown that no systematic offset between the synthetic isogravities and 
the observational point positions exists. This good consistency really confirms that the synthetic 
and the observational spectra are on the same 9 = 5040/T e j^ temperature scale and allows us to be 
confident on the methodology adopted for the comparisons. In general, for the dwarfs cooler than 
6000 K, the good agreement of the observational point positions with the synthetic isogravities in 
each metallicity bin convincingly indicates the correctness of the behavior of our defined synthetic 
indices versus surface gravity. Furthermore, the good consistencies in all [Fe/H] bins show that the 
dependency of the observational indices on metallicity could be well reproduced by our synthetic 
indices. However, a powerful evidence for agreement at very low metallicity ([Fe/H] < -1.0) and/or 
at effective temperature lower than 4200 K {6 > 1.2) is unfortunately hampered by few suitable 
data available from the test sample of the INDO-U.S. library. 

As for giants, we similarly plot the spectral index Cal, Ca2, Ca3, and CaT versus 6 (6 = 
5040/T e ff) in the eight different metallicity bins overplotted by SSA and NSSA isogravities of the 
calibrated indices respectively in Figure 7, 8, 9, and 10. In these figures, there seems to be a 
slightly tendency of the four synthetic indices to be overestimated in all metallicity bins. This less 
good consistency between the positions of the observational points and the synthetic isogravities 
most probably stems from the inaccuracy of stellar parameters for the INDO-U.S. giants because 
the consistency for giants is even much worse if we replace the improved stellar param eters from 



Wu et al.l (|201ll ) by the previous parameters from the INDO-U.S. library itself. Although IWu et aL 
'ias refined the stellar parameters for the INDO-U.S. library, the accuracy is still not enough 



for giants. Still, a sound comparison at effective temperature larger than 5000 K and at [Fe/H] < 
-1.0 and [Fe/H] = 0.5 is hampered by few suitable data from the test sample of the INDO-U.S. 
library. 

On the whole, our calibrated synthetic indices can give a good representation of the obser- 
vational indices, and such an representation is much better for dwarfs. Simultaneously, we are 
looking forward to a more precise version of the stellar parameters for the INDO-U.S. library, and 
particularly for giants. 



3.2. Test the calibrated synthetic indices on the SDSS stars 

In this section, we compare the calibrated synthetic Can indices with those measured on the 
SDSS-DR7 and SDSS-DR8 sample. A sample of 11,938 SDSS-DR7 stars (8309 dwarfs and 3629 
giants) with atmospheric parameters within our synthetic grid (3500 K ^ T e jf ^ 7500 K, 0.5 dex 
^ log g ^ 5.0 dex and -2.6 dex ^ [Fe/H] ^ 0.6 dex) is obtained from SDSS-DR7 spectroscopic 
database. Here, the "dwarfs" stands for those with surface gravities larger than 3.5 dex while stars 
with surface gravities smaller than 3.5 dex belong to the "giants". In addition, another sample of 
22,709 spectra of SDSS-DR8 stars (12,497 dwarfs and 10,212 giants) with parameters within the 
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Fig. 3. — Spectral index Cal of observed dwarfs stars vs. 9 = 5040/T e yy in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
black, green, blue, and red color represent log g = 3.5, 4.0, 4.5, and 5.0, respectively. 
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Fig. 4. — Spectral index Ca2 of observed dwarfs stars vs. 9 = 5040/T e yy in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
black, green, blue, and red color represent log g = 3.5, 4.0, 4.5, and 5.0, respectively. 
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Fig. 5. — Spectral index Ca3 of observed dwarfs stars vs. 9 = 5040/T e yy in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
black, green, blue, and red color represent log g = 3.5, 4.0, 4.5, and 5.0, respectively. 
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Fig. 6. — Spectral index CaT of observed dwarfs stars vs. = 5040/T e jy in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
black, green, blue, and red color represent log g = 3.5, 4.0, 4.5, and 5.0, respectively. 



-14- 



[Fe/H] = -2.50000 



[Fe/H] = -2.00000 




1.0 1.2 
5040/Teff 



[Fe/H] = - 1.50000 




1.0 1.2 
5040/Teff 



[Fe/H] = -0.500000 



4 

3 








o 2 

o 






♦ 1 


1 


i 















0.8 


1.0 


1 .2 1 .4 



5040/Teff 
[Fe/H]=0.200000 




1.0 1.2 
5040/Teff 




1.0 1.2 
5040/Teff 



[Fe/H] = - 1.00000 




1.0 1.2 
5040/Teff 

[Fe/H] = 0.00000 




1.0 1.2 
5040/Teff 

[Fe/H]=0.500000 




1.0 1.2 
5040/Teff 



Fig. 7. — Spectral index Cal of observed giant stars vs. 9 = 5040/T e j^ in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
orange, purple, dark grey, green, blue, and red color represent log g = 0.5, 1.0, 1.5, 2.0, 2.5, and 
3.0, respectively. 
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Fig. 8. — Spectral index Ca2 of observed giant stars vs. 8 = 5040/T e j^ in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
orange, purple, dark grey, green, blue, and red color represent log g = 0.5, 1.0, 1.5, 2.0, 2.5, and 
3.0, respectively. 
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Fig. 9. — Spectral index Ca3 of observed giant stars vs. 9 = 5040/T e j^ in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
orange, purple, dark grey, green, blue, and red color represent log g = 0.5, 1.0, 1.5, 2.0, 2.5, and 
3.0, respectively. 
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Fig. 10. — Spectral index CaT of observed giant stars vs. = 5040/T e yj in different metallicity 
groups. Synthetic SSA (solid lines) and NSSA (dashed lines) isogravities are overplotted. The 
orange, purple, dark grey, green, blue, and red color represent log g = 0.5, 1.0, 1.5, 2.0, 2.5, and 
3.0, respectively. 
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limits of our synthetic grid is also derived from the SDSS-DR8 spectroscopic database. 

Although the SNRs of our SDSS sample are all beyond 20, they are still much lower than those 
of the test sample from the INDO-U.S. library (Details about the effect of SNRs on the NIR Call 
indices are presented in Appendix Al). In addition to this, the number of the SDSS sample is quite 
large. So it is possible and better for us to show the tests by adopting the mean values instead 
of the individuals since it is always an effective way to diminish the effect of noises by adopting 
the mean value instead of the individuals under the circumstance that the number of individuals 
is quite large . Like what have been done in Section 3.1, we fill the eight different metallicity bins 
with the suitable spectra from our SDSS-DR7 sample, and in the next step, we group stars in each 
metallicity bin into several narrow bins by the effective temperature (9 = 5040/T e ^j) parameter at a 
constant step of A9 = 0.05. For now, we have gained numbers of sub-bins formed by the limitations 
of [Fe/H] and 9 together. Then, in each sub-bin where we firstly calculate the mean value of the 
effective temperatures and indices of all the stars, we could derive a couple of mean values (Meang 
and Meariindex) which absolutely represents a cursory generality of overall effective temperatures 
and indices in this sub-bin. Moreover, the mean values of the SSA and NSSA indices at every T e ff 
point from the synthetic grid are also computed in each of the eight metallicity bins. Then, we 
draw all couples of the means from all sub-bins in a Mearii nc i ex v.s. Meang panel, and for the sake 
of a clear comparison, the means of the SSA and NSSA indices are also overplotted as solid and 
dashed lines. Such comparisons for dwarfs and giants from the SDSS-DR7 sample are respectively 
shown in Figures 11 and 12 where every column presents the change of [Fe/H] from -2.5 dex to 0.5 
dex and every line exhibits the comparisons respectively for Cal, Ca2, Ca3, and CaT index. It is 
obvious that for all of the dwarf cases, the calibrated synthetic indices (solid and dashed lines) can 
well predict the trends of the indices from the SDSS-DR7 sample with 9 except for those in the 
bins centering on [Fe/H] = -2.5 and 0.5 dex where there are too few suitable data from the SDSS- 
DR7 sample. Detailedly, the agreement between synthetic indices and indices from the SDSS-DR7 
sample appears to be worse in the bin of [Fe/H]=-2.5±0.1 dex. As for the bin of [Fe/H]=0.5±0.1, 
there are no suitable data from the SDSS-DR7 sample. So we cannot be sure of anything conclusive 
for indices in these two metallicity bins. In order to solve this problem, we intend to remake the 
tests on the SDSS-DR8 sample which has more data in the bins of [Fe/H]=-2.5±0.1 and 0.5±0.1 
dex. Ultimately, the comparisons are shown in Figures 13 and 14 for dwarfs and giants. Evidently, 
the agreement between the calibrated synthetic and the observatianl SDSS-DR8 indices is also good 
for all the dwarf cases from the SDSS-DR8 sample but for those in the bins of [Fe/H]=-2.5±0.1 
and 0.5±0.1 dex where the SDSS-DR8 indices seem to be underestimated. While in the giant cases 
(Figures 12 and 14), we fail to see a good agreement between the synthetic and the SDSS indices 
because all of the SDSS indices fall below the theoretical lines which still most probably originate 
from the inaccuracy of the stellar parameters generated by the SSPP pipeline for SDSS stars. 

Conlusively, it is good to use SDSS dwarfs to check the validity of our synthetic indices; 
however, it is a little bad when SDSS giants are used. The reason is most likely to be the inaccuracies 
of the stellar parameters generated by the SSPP pipeline for the SDSS-DR7 and SDSS-DR8 giants. 
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So we hope for more precise stellar parameters for SDSS stars, especially for giants. Simultaneously, 
we are looking forward to the forthcoming LAMOST spectra which are expected to have more 
accurate stellar parameter estimations through a completely new and improved 2-D and 1-D stellar 
parametrization pipelines. 

4. Relations between Can indices and stellar atmospheric parameters 

4.1. Relations with surface gravity 

We investigate how dependent each of the indices (Cal, Ca2, Ca3, and CaT) is on the parameter 
of log g itself, and exhibit the investigation results in Figure 15 where the strength of an index 
is plotted versus log g for both the metal-poor ([Fe/H] < 0.0) and the metal-rich ([Fe/H] ^ 0.0). 
Apparently, for dwarfs (log g ^ 3.5; to the right region of the dashed line in Figure 15), there is 
little effect of log g on the four indices no matter whether the star is poor or rich in metallicity. 
Therefore, we could infer that, for dwarfs, the scatters of the indices in Figure 15 should be probably 
led by the effect of solely [Fe/H] or T e ff, or the combination of [Fe/H] and T e ff. Then, for giants 
(log g < 3.5; to the left region of the dashed line in Figure 15), indices strengths as a whole appear 
to be weakened by the increase of log g, which stands out for a general but negative influence of 
log g on index strength. Influence of such a kind looks more serious for the metal-rich than the 
metal-poor. Therefore, for giants also, the scatters of the indices at a same log g value should 
originate from the effect of solely [Fe/H] or T e ff, or the combination of [Fe/H] and T e ff. 

Furthermore, we find that indices strengths as a whole seem to be strengthened by the enhance- 
ment of [Fe/H] because the general level of the distributions of indices strengths of the metal-rich 
(the average line in the right pannel in Figure 15) is obviously higher than that of the metal-poor 
(the average line the left pannel in Figure 15), comparing the former with the latter. 

In short, it is revealed by Figure 15 that, for dwarfs, log g has a too little effect to be neglected 
on the indices strengths. Instead, [Fe/H] and T e ff play the dominant role; for giants, log g starts to 
affect indices strengths negatively which means that if log g increased, the indices strengths would 
decrease. As for [Fe/H], in general, the enhancement of metallicity can lead the indices strengths 
to be strengthened for both "dwarfs" and "giants" . 

4.2. Relations with metallicity 

We show the dependencies of Cal, Ca2, Ca3, and CaT index on the parameter of [Fe/H] in 
Figure 16 where the strength of each of the indices is plotted versus [Fe/H], respectively, for dwarfs 
and giants. Obviously, the four indices are being strengthened by the enhancement of [Fe/H] for 
both dwarfs and giants. However, for giants, in accord with Figure 16, the scatters of indices 
stengths at a fixed [Fe/H] appear more serious, and should probably stem from the effect of solely 
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Fig. 11. — Mean values of spectral index Cal, Ca2, Ca3 and CaT of SDSS-DR7 dwarfs v.s. 
6=5040/T e ff in different metallicity groups (columns). Mean synthetic SSA (solid lines) and NSSA 
(dashed lines) are overplotted for the clarity of agreement between our calibrated synthetic indices 
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Fig. 12. — Mean values of spectral index Cal, Ca2, Ca3 and CaT of SDSS-DR7 giants v.s. 
6=5040/T e ff in different metallicity groups (columns). Mean synthetic SSA (solid lines) and NSSA 
(dashed lines) are overplotted for the clarity of agreement between our calibrated synthetic indices 
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Fig. 13. — Mean values of spectral index Cal, Ca2, Ca3 and CaT of SDSS-DR8 dwarfs v.s. 
6=5040/T e ff in different metallicity groups (columns). Mean synthetic SSA (solid lines) and NSSA 
(dashed lines) are overplotted for the clarity of agreement between our calibrated synthetic indices 
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Fig. 14. — Mean values of spectral index Cal, Ca2, Ca3 and CaT of SDSS-DR8 giants v.s. 
6=5040/T e ff in different metallicity groups (columns). Mean synthetic SSA (solid lines) and NSSA 
(dashed lines) are overplotted for the clarity of agreement between our calibrated synthetic indices 
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Fig. 15. — Dependence of the Cal, Ca2, Ca3, and CaT index on the log g parameter. The calibrated 
LSS SSA indices are represented by the black points and the smoothed line represents the mean 
value which could give the general direction of the changes. The dashed line is a divide for dwarfs 
and giants. 
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log g or T e jf, or the combination of the two parameters. 

For dwarfs, considering the conlusion of log g having little effect on the four indices strengths 
from Figure 15, and then combining with Figure 16, we could easily infer that it is [Fe/H] and T e ff 
which affect the indices strengths together while log g plays an ignorable role. 

Then, for giants, [Fe/H], T e jj, and log g influence the indices together in strength. To be more 
detailed, for metal-rich giants, log g has a more serious effect on the indices strength than for the 
metal-poor giants, convinced by the much larger scatters of metal-rich giant indices stengths than 
those of the metal-poor giants in Figure 16. 

4.3. Relation with effective temperature 

We do try harder to find out how the parameter of T e ff affects the NIR Can indices via 
investigating how the indices stengths of SSA spectra change when T e jj changes. First, we divide 
the SSA spectra into four groups ( metal-rich dwarfs, metal-poor dwarfs, metal-rich giants, and 
metal-poor giants), and then measure the four indices on all groups. For each of the four indices, 
we calculate the average of index strengths over 9 (9 = 5040/T e jj) for each of the four groups, 
which are all depicted in Figure 17. Obviously, T e fj indeed has an effect on all indices, and only 
for the metal-rich giants presented by lines crossed with the 'plus' symbols in Figure 17, the effect 
appears obviously. 

To be summarized, the effect of the parameter of 9 [9 = 5040/T e jj) on Can indices is obvious 
for the metal-rich giants. For the sake of a clear but quanlitative illustration, the impact of the 
parameters of log g and [Fe/H] on the four indices is generalized in Table 2 where the 'little', 'much', 
or 'more' means that the parameter has little , much, or more effect on the four indices. The sign 
of '+' or '-' respectively stands for a positive or negative trend of the effect of the parameter on 
the four indices. 

4.4. Feature strength as a function of stellar atmospheric paramters 

In this section, we regress the relations of the index strength of CaT with stellar atmospheric 
parameters for synthetic stars, and then, test the derived relations on the observational stars 

Table 2: The qualitative effect of stellar parameters on Can indices. 



Para. 


Dwarfs 


Giants 




Metal-rich Metal-poor 


Metal-rich Metal-poor 


log g 


little little 


more (-) much (-) 


[Fe/H] 


much (+) much (+) 


more (+) more (+) 
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Fig. 16. — Dependence of the Cal, Ca2, Ca3, and CaT index on the [Fe/H] parameter. The 
calibrated LSS SSA indices are represented by the black points and the smooth line represents the 
mean value which could give the general trend of variations. The vertical dashed line is a divide of 
metal-poor and metal-rich stars. 
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respectively from INDO-U.S. library, SDSS-DR7, and SDSS-DR8 spectroscopic database which 
have already been used in the previous part of this paper. The regressions are made by the 
REGRESS program in IDL software. Suggested by the conclusions from Figures 15 ~ 17, for 
dwarfs, we expect to express the index strength of CaT by a function of metallicity ([Fe/H]) and 
effective temperature (9), and for giants, by a function of all the metallicity ([Fe/H]), surface gravity 
(log g), and effective temperature (9). The functions are below: 

For dwarfs (log g ^ 3.5) 

EW(CaT) = 3.648(±0.134)[Fe/H] - l.O81(±O.157)0 - 2.235(±O.119)[Fe/H]0 + 7.460 
For metal-poor giants (log g < 3.5 and [Fe/H] < 0) 

EW(CaT)=4.550(±0.211)[Fe/H] - 2.920(±0.176)log g - 6.765(±0.392)fl + 1.769(±0.147)log g 
9 - 1.966(±O.178)[Fe/H]0 - 0.354(±0.037)log g [Fe/H] + 17.450 

For metal-rich giants (log g < 3.5 and [Fe/H] ^ 0) 

EW(CaT)=11.955(±1.467)[Fe/H] - 5.491(±0.341)log g - 12.227(±0.646)# + 3.603(±0.297)log 
g 9 - 5.3O4(±1.234)[Fe/H]0 - 1.445(±0.258)log g [Fe/H] + 25.212 

The errors in the brackets are the 1-sigma uncertainty estimates for the corresponding co- 
efficients. The coefficients and the corresponding 1-sigma errors are directly returned from the 
REGRESS program. 

Using the functions above, we predict the strengths of CaT index on dwarfs, metal-poor gi- 
ants, and metal-rich giants respectively from the previously used SSA, INDO-U.S., SDSS-DR7, 
and SDSS-DR8 spectral sample. Then we compare predictions with the index strengths directly 
from measurements. For dwarfs, metal-poor giants, and metal-rich giants, the comparisons are 
respectively shown in Figures 18, 19, and 20 where the 'sigma' in every panel presents the standard 
deviation of differences between predictions and measurments. 

Previously in Section 3, the agreement between the observed Call indices and our synthetic 
ones is not so good for giants as that for dwarfs, which inevitably causes the different precisions of 
the functions above. From Figure 18, the analytical relation of the strength of the CaT index with 
the stellar parameters for dwarfs is good for both the synthetic SSA dwarfs and the INDO-U.S. 
dwarfs with small 1-sigma errors. As for SDSS stars, the 1-sigma error turns a little larger for DR7 
and DR8 stars because we could not obtain the more accurate parameters for SDSS stars so that 
the parameters we have to use are from SSPP which are only estimations with a low accuracy. We 
should have a much better consistency between predictions and measurements for SDSS stars if 
more accurate parameters could be obtained. The distribution has a up-shift from the gradient 1 
line for SDSS-DR8 indices which might be caused by a systematic difference of the different pipeline 
for data reduction from that for SDSS-DR7 spectra. In Figure 19 for the metal-poor giants, the 
function is still good for the SSA spectra while it still works for the INDO-U.S., SDSS-DR7, and 
SDSS-DR8 spectra despite a small deviation when the strength of the CaT index is smaller than 
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5. In Figure 20 for the metal-rich giants, the function seems good for the SSA, INDO-U.S., and 
even SDSS-DR7 sample while it does not work well for the SDSS-DR8 ones with the strength of 
the CaT index smaller than 5. Generally speaking, our functions of estimating the strength of the 
CaT index with the stellar parameters works for dwarfs, metal-poor giants, and metal-rich giants 
of late-types with parameters (3500 < T eff < 7500 K, -2.5 < [Fe/H] < 0.5 and 0.5 < log g < 5.0) as 
accurate as possible. These functions are simple and approximately good since the 1-sigma errors 
can be accepted. 



5. Comparisons with the CaHK index 



The CaHK index originating from the Can Fraunhofer H and K features (Ca H(A 3936) and 
K (A 3934)) can also b e use d to measure calc i um ab undance in galaxies. This index has been defined 



m 



Serven et al.l (|2005l ) and iFranchini et al.l (|2010l ). respectively. The Ca H and K doublet features 



are strong in the integrated-light spectra of cool stars and elliptical galaxies largely composed of 
low-mass and cooler stars. 

Comparing with Can triplet, Ca H and K give much stronger percentage response to increased 
Ca abundance, and the response spreads over a much broader wavel e ngth region which could 
make measurements even easier (details in Figure 1 in IWorthev et al.l (|201ll )). Furthermore, it 
i s noted that the CaHK index is a factor of three more sensitive than the Can triplet indices 
(jServen et al.ll2005l ). Unlike the Can triplet indices, the CaHK index is free of the contamination by 
the Paschen lines and the NLTE impact (details in Sections 2.1 and 2.2). However, the soundness 
of using the CaHK index as an ideal calcium measurement tool is hampered by the emission 
componen ts often present within the centers of th e H and K absorption lines for stars of type GO 
and later ( Schwarzschild et al. 1913 ). In specific, Wilson Sz Vainu Bappu (jl957l ) has furthermore 
found out the linear correlation of the logarithms of the emission-line widths with the spectroscopic 
absolute magnitudes. This Wilson-Bappu effect surely leads to weaker measurements of Ca H and 
K equivalent width but has no effects on measurements of Call triplet widths. Additionally, the 
Call triplet features appear around 8600 A in the NIR wavelength band where the absorbing effects 
of dust are small; instead, the Ca H and K doublet features occur around 3900 A in the very blue 
wavelength band where indices are more seriously subject to absorption by the interstellar medium. 



In 



Franchini et al.1 (|2010l ). it is shown that the CaHK index is sensiti ye to q-element enhanc e- 



ment which influences the CaHK index by strengthening it (Figures 1 in IFranchini et al.l (|2010l )). 
In our work, we have verified that the Can triplet indices are almost independent on a-element 
enhancement (Figure 1). For dwarfs, the behaviors of Call triplet indices with stellar parameters 
have a little better reproduced those of observational ones as can be deduced by the more compact 
consistency with the isograviti es in Figures 3 ~ 6 while the behaviors of the CaHK index are not so 



good as observed (Figure 6 in IFranchini et al.1 (J201Cj)). For giants, the Call triplet indices perform 
still a little better behaviors with stellar parameters (Figures 7~ 10) than the CaHK index (Figure 
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Fig. 17. — Indices of the Cal, Ca2, Ca3, and CaT versus 6 = 5040/T e jj. Lines crossed with 
symbols of star, triangle, plus, and square respectively represent the mean indices values of the 
metal-rich dwarfs ([Fe/H] ^ 0.0, log g ^ 3.5), the metal-poor dwarfs ([Fe/H] < 0.0, log g ^ 3.5), 
the metal-rich giants ([Fe/H] ^ 0.0, log g < 3.5), and the metal-poor giants ([Fe/H] < 0.0, log g < 
3.5). 




Fig. 18. — Predicted versus measured EW(CaT) for dwarfs. The predicted values have been 
calculated as functions of atmospheric parameters. The straight line has a gradient of 1. 
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Fig. 19. — Predicted versus measured EW(CaT) for metal-poor giants. The predicted values have 
been calculated as functions of atmospheric parameters. The straight line has a gradient of 1. 
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7 in iFranchini et al.l (|2010l )). Moreover, the light of the Ca H and K features is mostly contributed 
by main-sequence and warmer giants, in contrast to the light of the Call triplet featu res, which is 
dominated by red giant branch tip and asymptotic giant branch (jWorthev et al,ll201ll ). 

Briefly, the Ca H and K doublet are important features in the blue wavelength band, and so 
are the Can triplet in the NIR wavelength band. Both of them can be used together in spectra 
where both of the CaHK and the Call index are prominent in order to get a more convincing and 
consistent conclusion; they can be both useful tools to measure calcium abundance in galaxies under 
different circumstances and for different aims; however, once any of the two indices is adopted, it 
should be better to take its advantages and disadvantages into considerations. 



6. Summary 

In this work, we define a new set of NIR Can triplet indices and a new index CaT as the sum 
strength of the triplet lines. Then we measure every index of the four on the synthetic spectra, 
with the same spectral resolution (R = 1800) and wavelength scale (approximately 2.0) as SDSS 
and the forthcoming LAMOST spectra, which are synthesized by SPECTRUM, a spectral synthetic 
program, using the ATLAS9 model grid of which the effective temperature (T e fj) parameter ranges 
from 3500 to 7500 K with a step of 250 K, and the surface gravity (log g) covers from 0.5 to 5.0 dex 
with an interval of 0.5 dex. As for metallicity, two assumptions are adopted to model the chemical 
composition, Solar Scaled Abundances (SSA) with [a/Fe] = and Non-Solar Scaled Abundances 
(NSSA) with enhanced a element ([a/Fe] = 0.4 dex). In the SSA model, 8 values are adopted for 
metallicity ([Fe/H] = -2.5, -2.0, -1.5, -1.0, -0.5, 0.0, 0.2, 0.5) while in the NSSA model, a poorer 
value of [Fe/H] = -4.0 dex is supplemented. After measurement, we compare the Can triplet 
indices from the SSA model with those from the NSSA one at the same stellar parameters. Such 
a comparison indicates that the Can triplet indices are all quasi-independent on a element. Later 
on, we calibrate our synthetic Call triplet indices with the observed stars of the INDO-U.S. library. 
The most important part of this work is the comparisons of our calibrated synthetic indices with 




Fig. 20. — Predicted versus measured EW(CaT) for metal-rich giants. The predicted values have 
been calculated as functions of atmospheric parameters. The straight line has a gradient of 1. 
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the observational ones which can be summarized as follows. 

1. We make this comparison between our calibrated synthetic Can indices and the observa- 
tional ones of INDO-U.S. and SDSS DR7 and DR8 stars with the parameters within the synthetic 
grid. Comparison is respectively made for dwarf and giant stars. 

2. For dwarfs, the behaviours of our calibrated synthetic indices versus effective temperature, 
surface gravity and metallicity well reproduce that of the observational indices of INDO-U.S., SDSS 
DR7 and DR8 stars. It is worthwhile to stress that this good agreement is due to the fact that 
both theoretical and observational indices are on the same temperature scale. 

3. For giants, the agreement between the synthetic indices and the observational indices is 
not so good as that for dwarfs. This discrepancy might be the result of our inaccurate calibration 
caused by a little larger uncertainty of stellar parameters for giants. 

Therefore, we present a synthetic library of the NIR Can triplet indices and a summed index 
CaT which is quite appropriate and more trustful for fully exploiting the near-infrared content of 
the SDSS and the forthcoming LAMOST stellar spectroscopic database which represents and will 
represent the most extended collection of stellar spectra up to now. 

The relations of the four Call indices with stellar atmospheric parameters (log g, [Fe/H] and 
T e ff) have also been investigated qualitatively in the second half of this paper. Firstly, surface 
gravity (log g) has little effect on Can indices of dwarfs, but much effect of a negetive trend on 
those of metal-poor giants and even more on those of metal-rich giants. Secondly, metallicity 
([Fe/H]) affects much on the Can indices of dwarfs and more on those of giants by a positive trend. 
Speaking of effective temperature (T e ff), 9 (6 = 5040 /T e ff) certainly has an effect on the Call 
indices of both dwarfs and giants but by an irregular trend. Ultimately, we derive an function to 
approximately estimate the strength of the CaT index. The function for dwarfs well expresses the 
relation of the CaT index with the T e fj and [Fe/H] parameters and is tested on the SSA, INDO- 
U.S., SDSS-DR7, and SDSS-DR8 dwarfs. The predictions by the function are generally consistent 
with the measurements within the sigma errors which we could be accepted. All of the functions are 
applying for stars with more accurate stellar parameter estimates. Absolutely, the gain of the more 
accurate parameters for stars particularity for giants from the SDSS database is the key to improve 
our test of Call indices on the SDSS giants. Therefore, a new method to improve the precision of 
stellar parameters for the SDSS stars is expected. Simultaneously, we are looking forward to a large 
number of LAMOST stellar spectra with refined parameters to check the validity of our synthetic 
library of the Can indices on LAMOST spectra. 

In addition, we compare the differences between the Can triplet and the Ca H and K doublet 
indices in several major aspects. 

Finally, as a supplement of this paper, we investigate the effect of SNR to the new set of Call 
indices in the Appendix, which shows that the impact of SNR on the indices is so weak that we 
can ignore the effect when SNRs are beyond 20. 
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A. Sensitivity of the Can triplet indices to SNR 

The database on which we define our NIR Can triplet indices is made up of the synthetic 
spectra which have no noises. In this Appendix, we discuss the sensitivity of Cal, Ca2, Ca3, and 
CaT index to the spectral SNR. First we regard the SNR between 8620 and 8630 A of a synthetic 
spectrum as the SNR of the whole synthetic spectrum because the continuum is fairly smooth and 
few features exist in that wavelength region where the flux level of the continuum is quite close to 
one since our synthetic spectra are already normalized. In a clean theoretical spectrum, the flux 
of a pixel is the pure signal of the pixel. So the strength of the signal of the smooth continuum 
in a normalized theoretical spectrum is undoubtfully approximate to one. Then we can derive 
any SNR value we expect by adding random normally distributed Gaussian noises with standard 
deviations equal to the corresponding values to the whole synthetic spectrum. For example, if 
we add a random Gaussian noise with the mean and standard deviation respectively equal to 1 
and 0.2, we can obtain a noise-added synthetic spectrum with SNR = 1/0.2 =5. If the standard 
deviation is chosen to be 0.05, the SNR goes to be 20 (SNR = 1/0.05). In this test, we derive 23 
noise-added synthetic spectrum with different SNR (SNR = 1, 3, 5, 10, 15, 20, 25, 30, 95, 100) 
by adding normally-distributed Gaussian noises with the distinguished standard deviation values 
to one synthetic spectrum. We compute all of the four indices on each of the 23 different-SNR 
synthetic spectrum (Ca*) and those on the original synthetic spectrum with no noise added to 
(Ca). We do the same for every synthetic spectrum of our LSS SSA. Finally, Ca* and Ca are 
compared for all the SNRs. The comparison shows that the four indices are insensitive to SNR 
if SNRs are larger than 25, a little sensitive if SNRs are between 15 and 25 and very sensitive 
if SNRs fall below 15. However, CaT index is the sum of the Can triplet indices which means 
noises are summed simultaneously. So it should be more sensitive to SNRs. It seems that CaT 
index is affected more by SNR when SNRs are below 40. For NSSA, we get the same conclusions. 
Therefore, the lower the SNR is, and the less trust the four indices deserve. Finally, Figures Al 
~ A4 are presented to respectively show the changes of sensitivities of Cal, Ca2, Ca3 and CaT to 
SNRs. In every figure, we just show comparisons when SNRs are 15, 20, 25, 30, 40, 60, 80 and 100. 

This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 21. — Sensitivity of Cal index to SNRs. Cal* represents indices from synthetic spectra with 
no noise while Cal means indices from noise-added synthetic spectra. The left pannel shows Cal* 
versus Cal. The middle pannel directly shows the differences between Cal* and Cal, and 'dev' 
written down in this pannel is the standard deviation of the scatter. The black line is the zero-level 
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Fig. 22. — Sensitivity of Ca2 index to SNRs. Ca2* represents indices from synthetic spectra with 
no noise while Ca2 means indices from noise-added synthetic spectra. The left pannel shows Ca2* 
versus Ca2. The middle pannel directly shows the differences between Ca2* and Ca2, and 'dev' 
written down in this pannel is the standard deviation of the scatter. The black line is the zero-level 
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Fig. 23. — Sensitivity of Ca3 index to SNRs. Cal* represents indices from synthetic spectra with 
no noise while Ca3 means indices from noise-added synthetic spectra. The left pannel shows Ca3* 
versus Ca3. The middle pannel directly shows the differences between Ca3* and Ca3, and 'dev' 
written down in this pannel is the standard deviation of the scatter. The black line is the zero-level 
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Fig. 24. — Sensitivity of CaT index to SNRs. CaT* represents indices from synthetic spectra with 
no noise while CaT means indices from noise-added synthetic spectra. The left pannel shows CaT* 
versus CaT. The middle pannel directly shows the differences between CaT* and CaT, and 'dev' 
written down in this pannel is the standard deviation of the scatter. The black line is the zero-level 



